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The Green- Schwarz anomaly-cancelling mechanism in string theories requires a Chern-Simons 
term in the Einstein-Hilbert action, which leads to an amplitude birefringence of spacetime for the 
propagation of gravitational waves. While the degree of birefringence may be intrinsically small, 
its effects on a gravitational wave will accumulate as the wave propagates. The proposed Laser 
Interferometer Space Antenna (LISA) will be sensitive enough to observe the gravitational waves 
from sources at cosmological distances great enough that interesting bounds on the Chern-Simons 
may be found. Here we evaluate the effect of a Chern-Simons induced spacetime birefringence to 
the propagation of gravitational waves from such systems. We find that gravitational waves from in 
i coalescing binary black hole system are imprinted with a signature of Chern-Simons gravity. This 

, signature appears as a time-dependent change in the apparent orientation of the binary's orbital 

angular momentum with respect to the observer line-of-sight, with the change magnitude reflecting 
f\j ' the integrated history of the Chern-Simons coupling over the worldline of a radiation wavefront. 

, While spin-orbit coupling in the binary system will also lead to an evolution of the system's orbital 

■ angular momentum, the time dependence and other details of this real effect are different than the 

Q apparent effect produced by Chern-Simons birefringence, allowing the two effects to be separately 
indentified. In this way gravitational wave observations with LISA may thus provide our first and 

■ only opportunity to probe the quantum structure of spacetime over cosmological distances. 

7— I ' 

PACS numbers: 11.25.Wx, 95.55.Ym, 04.60.-m, 04.80.Cc 

cr 

bJQ, I. INTRODUCTION 

"Gravitational wave" is the name we give to a short-wavelength feature of the structure of spacetime, the arena 
within which all other phenomena play-out their roles. As such, the direct observation of gravitational waves offers 
an unprecedented opportunity to explore the environment that both enables and constrains the action of the broader 
^J. ' phenomena of nature. Here we describe an observational probe of quantum gravity made possible by the forthcoming 
, generation of space-based gravitational wave detectors. 

In most corners of the perturbative string theory moduli space (i.e. Type IIB, I, Heterotic) in four-dimensional 
compactifications, the addition of a parity-violating Chern-Simons term to the Einstein-Hilbert is required for math- 
ematical consistency [58j . Furthermore in the presence of the Ramond-Ramond scalar (D-instanton charge), the 
Chern-Simons term is induced in all string theories due to duality symmetries [l], 0] • The Chern-Simons term is not 
unique to string theory, but in fact it could also arise in Loop Quantum Gravity, where the coupling is not necessarily 
limited to small values. In the strong gravity sector of this framework, this term arises to ensure invariance under 
large gauge transformations of the Ashtekar connection variables [3[. This correction plays the same role in string 
theory that an analogous anomaly-canceling correction to the quantum-chromodynamics path integral plays in the 
Standard Model While the anomaly-cancelling field in the Standard Model case interacts with photons, leading 
to significant observational constraints on its coupling, the anomaly canceling term considered here only affects the 
gravitational sector of the theory and is mostly unconstrained by observation [!, H, 0] • 

The Chern-Simons correction to the Einstein-Hilbert action, Chern-Simons gravity [1] , has been recently studied in 
connection with tests both in the cosmological @> S Oil [HI and the weak-gravitational contexts @,B0]- ^ n the context 
of gravitational wave theory, Chern-Simons gravity leads to an amplitude birefringence of spacetime for gravitational 
wave propagation i- e -j a polarization dependent amplification/attenuation of wave amplitude with distance 

propagated. Observation of gravitational waves that have propagated over cosmolo gica l distances, such as will be 
possible with the Laser Interferometer Gravitational Wave Observatory (LISA) [LJ, |LJ, |l5[ , provide the opportunity to 
measure or bound the magnitude of the birefringence and, correspondingly, provide the first experimental constraints 



o 



- 1—^ 

X 



'Center for Fundamental Theory, Institute for Gravitational Physics and Geometry, Department of Physics 

t Center for Gravitational Wave Physics, Institute for Gravitational Physics and Geometry, Department of Physics, Department of As- 
tronomy and Astrophysics 

i Center for Gravitational Wave Physics, Institute for Gravitational Physics and Geometry, Department of Physics 



2 



on string theory models of gravity. 

Gravitational wave observations have long been recognized as a tool for testing our understanding of gravity (see [l6j] 
for a recent review) . Eardley and collaborator s (l7L LL8| first proposed a far-field test of all metric theories of gravity 
through gravitational wave observations. Finn [l9j|. and later Cutler and Lindblom [20| . proposed a means of realizing 
these measurements using a space-based detector in a circumsolar orbit observing solar oscillations in the far-zone field. 
Ryan 21;] argued that observations of the phase evolution of the gravitational waves emitted during the gravitational- 
wave driven inspiral of, e.g., a neutron star or stellar mass black hole into a supermassive black hole could be used to 
"map-out" the spacetime metric in the vicinity of the black hole horizon, testing the predictions of general relativity 
in the regime of strong-fields. There have been several proposals describing different ways in which grav itational wave 
observations could be used to place bounds on the graviton Compton- wavelength [12, [H, [H, Hf| [26ll27j , the existence 
of a scalar component to the gravitational interaction [2^, [2(| H3, HH , and the existence of other corrections to general 
relativity as manifest in some fundamental, dimensionful length scale [29l. [30|. The measurements we propose here 
are, we believe, the first example of a direct model-independent probe of string theory and quantum gravity with 
gravitational waves. 

In section |n] we review Chcrn-Simons modified gravity, focusing attention on the scale of the Chern-Simons term 
and its effect on the propagation of gravitational waves in a cosmological background. In section HTT1 we evaluate the 
observational consequences of the Chern-Simons term in the context of ground- and space-based gravitational wave 
detectors. In section [V] we summarize our conclusions and discuss avenues of future research. 

Conventions used in relativity work and conventions used in quantum field theory work are often at odds. We follow 
the relativity conventions Misner, Thorne and Wheeler [3l| in this work: in particular, 

• Our metric has signature — h + + ; 

• We label indices on spacetime tensors with greek characters and use latin indices to label indices on tensors 
defined on spacelike slices; 

• We use a semicolon in an index list to denote a covariant derivative (i.e., VyU becomes V V U KV ) and a comma 
to denote ordinary partial derivatives; 

• Except where explicitly noted we work in geometric units, wherein G = c = 1 for Newtonian gravitational 
constant G speed of light c. 

Note that in geometric units, units of mass and length are interchangeable (i.e., G/c 2 has units of (length) /(mass)). 
This is in contrast to Planck units (h — c — 1), where units of mass and units of inverse length are interchangeable 
(i.e., h/c has units of (mass) x (length)). 



II. CHERN-SIMONS MODIFIED GRAVITY 



A. Brief Review 



In this subsection we review the modification to classical General Relativity by the inclusion of a Chern-Simons 
term, based on [l,[l2]. All four dimensional compactifications of string theory lead, via the Green-Schwarz anomaly 
canceling mechanism, to the presence of a four dimensional gravitational Chern-Simons term Chern-Simons forms 
are formally defined for odd dimensions, with the 3-form of particular interest for gauge theories. By introducing 
an embedding coordinate, which may be dynamical, Jackiw and Yi Q described a Chern-Simons correction to the 
Einstein-Hubert action: 

S cs - ^Jd 4 x6R*R (2.1) 

where 9 is (a functional of) the embedding coordinate, 

R*R=^R aPlS e a ^W\ lv , (2.2) 

and e a P"<5 j s the Levi-Civita tensor density. The variation with respect to the metric of this contribution to the total 
action (which includes the Einstein- Hilbcrt action plus the action corresponding to any additional matter fields) yields 

5S CS = -jL f d 4 xV=g C af3 6g a P, (2.3) 
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where g is the determinant of the metric and C is the C-tensor [5 



£ja[3 



1 



e7<5(a 



(2.4) 



and the parenthesis in the superscript stand for symmetrization. The variation of SqSi the usual Einstein- Hilbert 
action, and the action of other matter fields leads to the equations of motion of Chern-Simons modified gravity: 



G a l3 + C a p — 8irT a p, 



(2.5) 



where G a b is the Einstein tensor (i.e., the trace-reversed Ricci tensor) and T t, is the stress-energy tensor of the matter 
fields. 

By construction the divergence of the Einstein tensor G a /3 vanishes. If 9 is treated as a fixed, external quantity 
then general covariance, which requires V ■ T = 0, leads to the constraint V ■ C = 0, which is shown in [8[ to be 
equivalent to R*R = 0. Alternatively, if 9 is a dynamical field, then variation of the action with respect to 9 will 
lead to the same constraint on R*R. Here we are interested in the propagation of gravitational waves in vacuum, 
where T = and the constraint V • C = is satisfied regardless of whether we view 9 as a dynamical field or a fixed, 
externally-specified quantity. 



B. Linearized Chern-Simons gravitational waves 

Focus attention on gravitational wave perturbations to a Friedmann-Robertson- Walker (FRW) cosmological back- 
ground in Chern-Simons gravity. Following [Tig , we can write the perturbed FRW line element as 

ds 2 = a 2 ( V ) [-dif + (Sij + h l3 )dx l dx 3 } (2.6) 

where rj is conformal time, x l are comoving spatial coordinates, Sij is the Euclidean metric, and hii is the metric 
perturbation, which — for gravitational wave solutions — we can take to be transverse and traceless [60| . Introducing 
this perturbation into the field equations [Eq. (|2.5D ] leads to 

U a h?i = -^ Pjk W - 2W) d p h' ki + 9' d p D g h kl } , (2.7) 
where we have introduced the notation 

' = d v (2.8) 

□ fl = d 2 n - didrf* + 2Hd v (2.9) 

H = a! ja. (2.10) 

Conformal time r\ is related to proper time measured by an observer at rest with respect to the cosmological fluid via 

dt = a(r))dri; (2.11) 

correspondingly, the conformal Hubble function 7i is related to the Hubble function H measured by an observer at 
rest with respect to slices of homogeneity via 

H=- = -H, (2.12) 
a a 

where we have use dots to stand for derivatives with respect to cosmic time t. 

Focus attention on plane-wave solutions to the wave equation [Eq. (|2.7[) ]. With the ansatz 

hi m (v,X l ) = ^e-^-^ k ], (2.13) 

where the amplitude Ai m , the unit vector in the direction of wave propagation nu and the conformal wavenumbcr 
k > are all constant, we find that 4>, k and Aij must satisfy 

VAij = -a- 2 e pjk n p A ki [(9" - 2H9') {<p' ~iH)n + i9'nV] (2.14) 
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where 

V= -i(f)"- {4>'f - H' - H 2 + k 2 . (2.15) 

Since the Chern-Simons correction breaks parity, it is convenient to resolve Aij into definite parity states, corre- 
sponding to radiation amplitude in the right- and left-handed polarizations efj and efj : 

Aij = Anefj+ALefj (2.16a) 



where 



1 
1 

3 



u - -=(4+ ie «) ( 2 - 16b ) 

- ^(4l - ie kl) ( 2 - 16c ) 



and eti x are the usual linear polarization tensors [3l| • It is straightforward to show that 



fe ^ L =*A R , L (e^) R ' L , (2.17a) 



riie ■> e 
where 

Ar = +1 (2.17b) 

A L = -1. (2.17c) 

With this substitution equation (|2.14[) becomes two decoupled equations, one for right-hand polarized waves and one 
for left-hand polarized waves: 



W + (0r,l) + W + U - = (i_A R , L «0'/a2) (2 ' 18) 



The terms on the right-hand side of equation (|2. 18[) are the Chern-Simons corrections to gravitational plane-wave 
propagation in a FRW spacetime. To understand the relative scale of these terms, we rewrite the equation in terms 
of the ratio <j>' / k: 

— +i(l-jT ~SA-y) = — (y-i<yT) (2.19a) 

where 

y = - (2.i9b) 

7 =— and T = ^ (2.19c) 
k no 

5=^- and A = ^- (2.19d) 

9" 9" 

e = 4 and £ = (2.19e) 

( = ^ and (2.19f) 

% a c 

and a subscript indicates the present-day value of the functions 9' , 0", 7i, W and a. 
If we assume that 9 and 7i evolve on cosmological timescales (i.e., /' ~ Hf) then 

e 2 „ (7C) 2 <<; 7 2 „ |j| (2 20) 

Treating the terms in e and 7^ as perturbations, write the solution to equation (|2.19ap as 

y = yo + «/o,i + iCvifl + ■■■ (2.21) 



5 



where j/o = <^o/ K * s the solution to the unperturbed equation [i.e., the dispersion relation in an FRW cosmology, given 
by equation (|2.19a[) with vanishing right-hand side]. The first corrections j/o.i an d 2/1,0 owing to the Chcrn-Simons 
terms satisfy 

2/o,i ^ 2iK?/oyo,i = ^R,hnEy (2.22a) 
y' l0 - 2wy yi,o = -2\ R xkTZij . (2.22b) 

Requiring that the perturbation vanish at some initial (conformal) time rji the perturbations yo,i and 2/1,0 satisfy 

Ito.ifa) = A R , L y[-B](r?) (2.23a) 
yi,o(v) = -2A R>L 3>[rZ](rj) (2.23b) 



where 



y[ 5 ](r?) = /tc 2 ** ^ f 1 dxe- 2lMx) y a (x)g(x). (2.23c) 



Finally, the Chern-Simons correction to the accumulated phase as the plane-wave propagates from 77.; to 77 is 

^r,l = kAr,l /'d?7 {ey[^]( 7? )-2 7 Cy[rZ](7 7 )} (2.24) 



When 7 < 1, i.e., fco is very much greater than the Hubble constant Hq, the rescaled frequency \yo\ ~ 1. In this 
limit we can use integration by parts to find an asymptotic expansion for y[g]: 



y[g](ri) 



LA 

2in J \yo dr] 



' x \n+l 



0[- (2.25) 



In the next section we explore the observational consequences of gravitational wave propagation in Chern-Simons 
gravity. 

III. OBSERVATIONAL CONSEQUENCES 

A. Birefringence in a Matter-Dominated Cosmology 

Current and proposed ground-based gravitational wave detectors are sensitive to gravitational waves in the 10 Hz 
- 1 KHz band [32], [H, [34|, HE]- Detectable sources in this band are expected to have redshifts z < 1. Space-based 
gravitational wave detectors like LISA [l3[ will be sensitive to gravitational waves in the 0.1 - 100 mHz band and, 
in this band, be sensitive enough to observe the gravitational waves from the inspiral of several ~ 10 6 M Q black hole 
binary systems at z < 30: i.e., anywhere in the universe they are expected [36l. [37j. For sources in the band of these 
detectors 



where 

Ho 

100 km s _1 Mpc 



hwo — 7^7-; _i^ T _i ■ (3-2) 



Additionally, for redshitfs z < 30 the universe is well described by a matter-dominated FRW cosmological model. 
In this section we evaluate the effect that the Chern-Simons corrections described above have on propagation of 
gravitational plane-waves through a matter-dominated FRW model. 
In a matter-dominated FRW model the scale factor a(t) satisfies (HI 

^W = -L- (3.3a) 
a Q 1 + z 
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where, by convention, rj — 1 at the present epoch. In this model and with this convention 

2 



H=- = 2yTT^ 





j 


7 = 










2 


5= - 


^2 


e = M 
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O'a 


_ H Q k 
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Additionally, 



E 
Z 



and 
and 
and 
and 
and 



7^ 



H = 2, 



r = ir 1 = VT+z, 

A = if 2 = 1 + z, 
F = 1 9 1 

V 4 O'o' 



a = 



(3.3b) 

(3.3c) 
(3.3d) 
(3.3e) 

(3.3f) 
(3-3g) 

(3.3h) 
(3.3i) 
(3-3j) 



With the 7 and T parameters for a matter-dominated FRW cosmological model, the unperturbed equation for yo 
has solutions of the form 

i (l + Cnrj — K 2 ry 2 ) cos (ki]) — (C — nr\ — CK 2 -q 2 ) sin (^77) ^ 

k?7 (1 + CK77) cos (kt;) — (C — ktj) sin (k7/) 

where C is a constant of the integration. In the limit of large k the evolution of n should decouple from the universal 
expansion; thus, we are led to choose C — ±i, which eliminates the oscillatory terms in our general solution for yo: 

±K?if - i 



^° KT] (1 + n 2 r] 2 ) 

Consistent with our ansatz [cf. Eq. (|2.13p ] we choose the solution with positive 5i(i/o) : i.e., C 
equation for the phase 4>' = nyo we find 



&Mv) = & (1)- to fa) 



K (1 — tj) — arctan 



«(!->?) 
1 + K 2 rj 



- - In 
2 



1 + K 2 7] 2 

rf (1 + k 2 ) 



2 1 



7 + 2 



arctan 



2 7 (vm - 1) 

4 + 7 2 VTTI 



I In 
2 



2 

7 z z 



4 + 7 2 



(3.5) 
-i. Solving this 

(3.6a) 
(3.6b) 

(3.6c) 



In the absence of the Chern-Simons correction an observer at rest with respect to slices of (cosmological) homogeneity 
will observe a passing gravitational plane- wave to undergo a change in phase A0o( J 7) between cosmological time 77 
and the present epoch. 

With yo and equations (|2.23p we can evaluate the Chern-Simons contribution to the phase change owing to prop- 
agation from cosmological time 770- Making use of the asymptotic expansions for yo,i and yi.o [Eq. (|2.25p ] we find 



3/0,l(R,L) 
Z/l,0(R,L) 



2 

-«A RL 



E(T)) - £(^ o ) e -2z[A0 o ( t)o )-A<Mr ; )]j + {"f) 

T{rj)Z{ri) - r(?7o)Z(r7o)e^ 2l[A0o(?,o) - A0o(, ' )1 l +G( 7 ), 



which may be integrated to find A0x(r,l): 

A01(R, L ) 



zAr ' l 2 
Ac 



'(r,) 



rf O ' 



2 7 C 



1 



(3.7a) 
(3.7b) 

(3.8a) 
(3.8b) 
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It is convenient to rewrite £ as a function of z: 

£(z) = aA(z)+(3B{z), (3.8c) 

where 

m (3 - 8d) 

3. (dg/dz) 

7 ^ 2 2(d 2 0/dz 2 ) o + 3(^/dz) o ' 1 ; 
e(d 2 9/dz 2 ) Q 

13 ~ 2(d 2 d/dz 2 ) Q + 3(d0/dz) o ' (3,8g) 

and the subscript zero denotes present-day values of the subscripted quantities. The leading order Chern-Simons 
correction to the accumulated phase is thus pure imaginary, corresponding to an attenuation of one circular polar- 
ization state and an equal amplification of the other. The attenuation/amplification is linearly dependent on the 
wavenumber. The function £(z) may be thought of as a "form-factor" that probes the past history of the coupling 6. 



B. Binary Inspiral at Cosmological Distances 



The proposed LISA gravitational wave detector is capable of observing coalescing binary black hole systems at 
cosmological distances: for example, the gravitational waves associated with a pair of 1O 6 M0 black holes will be 
observable at redshifts z approaching 30. Over the year leading up to the merger of two such black holes the binary's 
period will decrease by two orders of magnitude, leading to a corresponding decrease in the radiation wavelength and 
increase in the magnitude of the Chern-Simons correction. The time-dependent relationship between the radiation 
amplitude in the two polarization states thus carries with it the signature of Chern-Simons gravity and can be used 
to characterize the functional 8 that describes the Chern-Simons correction to classical General Relativity. 

To calculate the signature left by the Chern-Simons correction on the gravitational waves from a coalescing binary 
system at redshift z, we begin with the radiation near the source. Treating, as before, the Chern-Simons correction 
as a perturbation, the quadrupole approximation to the radiation from the binary system in the neighborhood of the 
source is given by 



h = 3? h + e + + h x e x 

k = ^[i+x 2 ] [Mm/* 



2/3 



exp 



Mk{t)/2 



2/3 



exp 



i ($(£) - k{i)d 
-i ($(t) - jfc(t)d) 



where d is the proper distance to the source and 

<£(£) = -2 



k(t) = ^ 




M \ 256r-t 



(3.9a) 
(3.9b) 
(3.9c) 

(3.10) 
(3.11) 



The constants f and 5, which determine when coalescence occurs and the phase of the gravitational wave signal at 
some fiducial instant, are set by initial conditions. The quantities M and x are constants that depend on the binary 
system's component masses (mi, m-i) and orientation with respect to the observer: 



X 
M 



cosine-angle between the orbital angular 
momentum and the observer line-of-sight 



3/5 3/5 
TO]/ TO 2 



(mi + m 2 ) 



1/5 



= ( "chirp" mass) . 



(3.12a) 
(3.12b) 
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We "hat" all these quantities to remind us that, as expressed above, they are appropriate descriptions only in the 
neighborhood of the source where the Chern-Simons and cosmological corrections to the propagation of the waves 
may be neglected. 

To describe the radiation after it has propagated to the detector we first describe the near-source radiation in terms 
of circular polarization states 



IR.L 



V2 



hR e R + h L e L 
M ( Mk\ 2/i 



(1 + Ah : lx) 2 exp -i ( $(£) - k{t)d 



(3.13a) 
(3.13b) 



We are interested in the radiation incident on our detector today (z = 0, r\ = 1) from a source at redshift z. Matching 
the near-source description of the radiation [Eq. (|3.13|) ) to our ansatz [Eq. (|2.13j) ] we find the description of the 
radiation after propagating to the detector from a redshift z: 



3? 



fiR e R + h L e L 

2/3 



where 



*o(i) = 

k (t) = 

d L = 

M = 

A^o = 



lit) 



T-t 



5/8 



5M 
2/5 M 



\-s, 

3/8 



M V256T-t / 
aor](l + z) — (Luminosity distance to source) 

(l + z)M 
2(VTTz- l) 



-k (t) ~ arctan27(<) 



VT+z-1 



4 + 



In 



1 + 7 2 W 



4 + 7 2 (i) 



Hp 
k a (t)' 



(3.14a) 
(3.14b) 

(3.14c) 

(3.14d) 

(3.14e) 
(3.14f) 

(3.14g) 
(3.14h) 



and A^imH given by equations (|3.8p above. Here t is proper time as measured by a detector at rest with respect to 
the cosmological fluid at the present epoch (rj = 1), ko(t) is the instantaneous wavenumber of the wavefront passing the 
detector at observer time t, and T and i5 are, as before, constants of the integration. The correction A^o, which is the 
same for all polarizations, embodies O(ko/Ho) corrections to the wave phase owing to the wave propagation through 
the time-dependent cosmological background. The correction A^mn is of opposite character for the two polarization 
states and embodies the (first-order) corrections to wave propagation owing to the Chern-Simons corrections to the 
Einstein Field Equations. 

Focus attention on the argument of the exponential in equation (|3.14b[) . The term k (1 — rj) cancels the first term 
in equation Q3.14gp for A0o, leading to 



1R,L 



V2 



M fMk \ 2/3 



(1 + A fl ,LX) exp 



-'■(*o(t)-^(VT 



A(f> 



1(R,L 



(0 



(3.15) 



The observational effect of the Chern-Simons is readily identified by looking at the ratio of the polarization amplitudes 
h,R and h^: 



(1 + X) 

(i-x) 

1 + x 
l-x 



exp 



Ho 



(3.16) 
(3.17) 
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where £ is given by equations (|3 - 8[) and x may be interpreted as the apparent inclination cosine-angle. The effect of 
the Chern-Simons correction on gravitational wave propagation is to confound the identification between polarization 
amplitude ratios and binary orbit inclination cosine-angle. In the same way that we say that the curvature of spacetime 
"bends" light passing close to strongly gravitating body we may say that the effect of the Chern-Simons correction is to 
"rotate" the apparent inclination angle of the binary system's orbital angular momentum axis either toward or away 
from us. 

IV. DISCUSSION 

A. What can be measured? 

Over the course of a year-long observation the LISA spacecraft constellation will measure the radiation in both 
polarizations of an incident gravitational wave train associated with an inspiraling coalescing binary system. The 
relative amplitude of the two polarization will be determined by the orientation of the binary systems orbital plane 
to the observer line-of-sight and the form factor £(z). A non- vanishing £ leads to a time- varying apparent inclination 
angle that, by nature of its time dependence, can (in principle) be measured directly from the apparent inclination 
angle's time variation. 

Other properties of an inspiraling binary can lead to an evolution of the (apparent) inclination cos-angle. Spin-orbit 
coupling leads to a real precession of the binary's orbital plane and a corresponding time-dependence in the actual 
inclination cos-angle \. Referring to equation (|3.14b[) , it is apparent that for small \\\ ~ an incremental change /j in 
y w ill lead to changes in /ir^l that are indistinguishable from an increment in x associated with £. Following Vecchio 
[3« . eq. 27-31] we note that, at first non-vanishing post-Newtonian order, spin-orbit interactions in an inspiraling 
binary system lead to 

(§)_«#*<*>• («) 

This is a different dependence on ko than the O(fco) dependence associated with £. Thus, it remains in principle 
possible to distinguish the signature of Chern-Simons gravity in the signal from cosmologically distant coalescing 
binary black hole systems. The accuracy with which such a measurement can be made is the topic of the next 
subsection. 



B. How accurately can £ be measured? 

The most general astrophysical black hole binary system can be described by eleven independent parameters, which 
may be counted as two component masses; component spins and their orientation (six parameters); orbital eccentricity; 
orbital phase; and a a reference time when the phase, spins and eccentricity are measured. The gravitational wave 
signal in any particular polarization will depend on the description of the binary and six additional parameters that 
describe the binary's orientation with respect to the detector. These six additional parameters may be counted as 
orbital plane orientation (two angles); source location with respect to the detector (distance and two position angles); 
orbit orientation in orbital plane (one angle) [6lj . To these seventeen parameters we now add £, which describes the 
effect of propagation through the birefringent Chern-Simons spacetime, for a total of eighteen parameters that are 
required to describe the signal from a coalescing binary system. 

To-date, all analyses of expected parameter estimation errors have been made with under a set of approximations 
that focus attention on the measurement of component masses, source location (both distance and angular position), 
and the expected time of binary coalescence. Even the most sophisticated of these analysis i gno re all but the leading- 
order contribution to the gravitational wave signal amplitude at twice the orbital frequency [62| and assume that the 
orbital eccentricity is known to vanish. These approximations are quite appropriate for their purpose (estimation of 
component masses, source location and expected time of coalescence); however, by ignoring all but the leading order 
contribution to the signal magnitude they are inadequate starting points for exploring the accuracy with which £, 
which affects only the signal amplitude in the different polarizations, can be bounded [63]. Evaluating and presenting 
the errors associated with the measurement of £ via a full co-variance matrix analysis is thus a formidable enterprise, 
to be addressed in a future work. 

Nevertheless, through a series of plausible approximations it is possible to make a crude estimate of the accuracy 
with which £ can be determined. To begin, assume we have two gravitational wave detectors such that, via a linear 



10 



combination of observations made at each, we can can synthesize two other detectors with one exclusively sensitive 
to /ir and and one exclusively sensitive to hi,. Write the scalar detector response of each of these detectors as 



m-R,x{t) = exp [/^r,l (t) + #r,l(*)] 



(4.2) 



for real /^r,l and V'R.l- Next, note that the parameters that describe a coalescing binary system can be divided into 
two groups: those that principally affect only the signal amplitude (i.e., /i(i)) and those that affect only or principally 
the real part of the signal phase (i.e., ip(t)). The first group includes distance, source orientation with respect to the 
observer line-of-sight, and £. The second group includes the orbital phase, sky location (through its affect on the 
Doppler correction to the signal phase as the detector orbits about the sun), the instantaneous binary period at some 
fiducial moment, and the parameters associated with spin and orbital angular momentum [64j |. If we approximate each 
detectors' noise as white with one-sided noise power spectral density So then the elements of the inverse covariance 
matrix T — the so-called Fisher matrix — are given by [H, H3| 



r — 



E 

fc=R,L 



2 

So" 



\ dx l 



7? 



dt 



(4.3) 



where the integration is over the observation period (ti,tf) and the x % are the parameters that characterize the incident 
gravitational wave, which we have divided into two groups. Matrix elements T^- where x l and x J belong to different 
groups will be much smaller than elements where x 1 and x J belong to the same group. Setting the cross-group elements 
to zero we obtain an approximate T that is block diagonal, with one block corresponding to ITy with (x l , x J ) drawn 
from the first group and the other block corresponding to I\y with (x l , x J ) drawn from the second group. Estimation 
uncertainties of parameters in either group can now be determined independently of the parameters in the other 
group. 

Focus attention now on those parameters that affect only n(t), the signals amplitude evolution. The leading order 
dependence of the amplitude |/ir,l| on the binary systems parameters is given by 



Ir.l = |/ir,l| = (1 + Ar,lXo) 2 



k (t)M 



2/3 



exp Ar,l£ 



kp(t) 
Ho 



(4.4) 



where M is assumed known. Setting aside the antenna pattern factors associated with the projection of the signal 
onto the LISA detector (which depend only on the known source sky position and the LISA orbital ephemeris), 
assuming that there is no real precession in the binary system under observation (i.e., xo), and that k(t) is given by 
equation (|3.11|) the inverse of the covariance matrix — the so-called Fisher matrix, T — associated with the amplitude 
measurements is a symmetric 3x3 matrix with elements 



- T>xo 



r-D£ — — 



x XoXo 



r X0« 



r« = 



l 
i 

Sb 

i 

So 
1 

So 
1 

So 
1 

So 



(A 



At) 



dt 



6Xo 



xt) 



M 



1 + 0(0 



1 -Xo 
k (t) 



Al 



Hn 



(At — A 



4^4r 



— A 2 

f Xo 

dt ~ - 



4At 



(l + xoY (i-xoY 



f 2k (t) 
H 

! k 2 (t) 
HI 



A 2 
1 + Xo 



Al 



dt ~ -16xo (3 

64xq (1 + xl] 
MH 

dt ~ 32 (l + 

32 (l 



xl) 



M 

dL 



dL ) 



xl) 



M 

d L 



(A 



Al 



1 

dt 



dt 



Zxt) 



Xo_ 
32 (1 



6xg 
(MH ) 



MH 

xt) 



o(0 
+ o(0 
o(0 



M 

d L 



K + O(0 



(4.5a) 
(4.5b) 
(4.5c) 
(4.5d) 
(4.5e) 
(4.5f) 
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where 



V = \nd L 

( k {t)M ^ 



4/3 dt 5 ' fc » 



J 
K 



"'f ( k (t)M \'- " <h r, 



7 So 192 

7/3 ^ K ,./,„ 

J So 192 

! ' / fc ft).M \ 10/3 dt _ _5_ 
J So ~ 192 



/fc7WV 7/3 7W 2 dfc 

S 32 
/fc7WV 1/3 7W 2 dfc _ 2 1 / 3 M 5 



2 1 / 3 7W 5 



(kM) 



V 2 / 



5 n 



128 



(feM) 



-1/3 



2/3 



(4.5g) 
(4.5h) 

(4.5i) 

(4-5j) 



and we have taken advantage of the fact that for inspiraling compact binary systems in the quadrupole approximation 
k(t) is monotonic in t to re-express the integrals over the interval (tj, tt) as integrals over [fc(ij), k(tt)] = (fc m inj &max)- 
In the particular case of a binary seen plane-on (xo = 0), the (T>xo) and (2?£) blocks of T are diagonal, leading to 



1 



Kl 



XnXo 



Ap 2 IK - J 2 



(MH a ) 



Ap 2 IK - J 2 



where z/jj is the ensemble average co- variance 



Uij = {x l - x t ) (xJ -xj) = (r ).. 
and we have expressed the in terms of the ensemble average amplitude-squared signal-to-noise ratio p 2 



P 2 = -/ dt(A% + A 



(4.6a) 
(4.6b) 

(4.6c) 
(4.6d) 



Focus attention on a binary system of two black holes at redshift z, each with mass M = 10 6 M Q (l + z) \ Over 
the final year before coalescence the radiation wavelength 2ir jk observed at the detector will range from c (l0~ 4 Hz) 1 
to c (l0~ 2 Hz) . For such a system, 



'100 



1 + Z- y/TTz 

= 3.1 x 1(T 40 



IP" 40 Hz' 1 
So 



1/2 



( — T ) (i + z- vThk£) s 

VlO- 40 Hz-V V ' 



(4.7) 
(4.8) 



Observation of binary systems like these at z = 15 by LISA will be capable of placing a "1-sigma" upper bound on £ 
of order 1CP 19 . 



C. How large might £ be? 



To estimate £ [cf. Eq. (|3.8cp ] we must invoke a theoretical model for the functional 9[(f>(z)]. As described in the 
introduction, perturbative string theory requires a Chern-Simons correction to the Einstein-Hilbert action [651 ]. Here 
we describe a different mechanism, that can also lead to the presence of a Chern-Simons correction. Consider the 
back-reaction of a J\f = 1 supersymmetric Yang-Mills theory in a curved background (cf. [ill . Appendix A]) with 
action 

Scs = ~ikJ dix • Fi(s) (R * R) ' (4 - 9) 

where S is the glueball superfield and .Fi(S), which plays the role of 9 in equation (|2.ip . can be exactly evaluated by 
using perturbative matrix model technology developed in ;42j. Within this Yang-Mills framework, 9 is a functional 
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9[ip] of some pseudo-scalar field ip, the gravitational axion, that depends only on conformal time [ll 
9[ip] can be expressed as 



Nil, 



2tt Mpi 



The functional 



(4.10) 



where £pj = (fi,G/c 3 ) is the Planck length, Mpi = (^pi-s/Stt) -1 is the reduced Planck mass and A/" is a dimensionless 
constant. Through use of the low-energy effective 4-dimensional heterotic string action the constant M can be 
evaluated in terms of the ten-dimensional fundamental string energy scale M s and the dimensionless string coupling 

9s, 



TV : 



( Mpi 



1 



leading to 



9[tp\ 



1 ip M] 



pi 



16^M^M B M s ^g s - 

Assuming that tp, which has units of inverse length, evolves with the Hubble parameter H oc ry -3 we have 



A(z) = B(z) 



2 r 
TT 



(l + z) 11 / 2 -! 



and 



with (e — 7f) of the order of 



(l + z) 11/2 -l (e-2 7 C) 



e - 2 7 C : 



(l.8 + 3.5fr? O0 ) x 10 



-120 



/5b 



\ / 10 16 GeV \ 
AfJ I. M s J 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15a) 



The size of £ thus depends on the present value of the field ipo, the fundamental string energy scale M s and the string 
coupling (7 S , none of which are constrained by present-day theory. 

The lesson to draw from the discussion of this scenario is that the magnitude of any Chern-Simons correction depends 
strongly on the external theoretical framework that prescribes the functional 0[ip]. For non- vanishing string coupling 
in the perturbative string theory scenario the Chern-Simons correction seems undetectable owing to the Planck scale 
suppression of the decay constant of the universal gravitational axion field ip. However, this model and the associated 
expected scale of e — 27^ applies only to the perturbative sector of string theory and, in particular, when Ramond- 
Ramond charges are turned off. If present these additional degrees of freedom do couple and source the Chern-Simons 
correction, leading to a larger decay constant (e.g, D3 branes always excite the Chern-Simons interaction in four 
dimensions). In a recent work, Svrcek (43| noted that, due to non-perturbative gravitational instanton corrections, 
the Chern-Simons coupling in the non-perturbative sector is currently incalculable. Even within the perturbative 
framework there are theoretical frameworks where £ co uld become significant: e.g., if the string coupling g s vanishes 
at late times IH, 0, H3, S, 0, HH [El, [H, [HI] . Therefore, within the full string theory framework, a larger 
coupling, which would push the stringy Chern-Simons correction into the observational window, is not excluded and 
bounding it places a constraint on string theory motivated corrections to classical general relativity. 

Setting aside string theory for a moment, other theoretical frameworks lead to a Chern-Simons correction to the 
Einstein-Hilbert action. Jackiw and Pi [8( showed that the embedding of the 3-dimensional Chern-Simons action into 4- 
dimcnsional General Relativity leads to Chern-Simons gravity as described by equation (|2.ip . Quantum mechanically, 
the Chern-Simons correction encodes information about the Immirzi parameter in loop quantum gravity and is 
also related to torsion. In this case, the Chern-Simons term is required to ensure invariance under large gauge 
transformations [3j. In both these classical and quantum mechanical scenarios the coupling of the Chern-Simons 
term to the Einstein-Hilbert action is theoretically unconstrained and is consistent with a coupling of order unity. 
Observations of gravitational waves from inspiraling binary black hole systems can thus be used to probe the strong 
gravity sector of quantum gravity. 
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V. CONCLUSIONS 

Chern-Simons type corrections to the Einstein-Hilbert action are strongly motivated by string theory, loop quantum 
gravity, and other scenarios. In all cases these corrections lead to an amplitude birefringence for gravitational wave 
propagating through space time. We have evaluated the correction to the gravitational wave amplitude for waves 
propagating over cosmological distances in a matter-dominated Friedmann-Robertson- Walker cosmology. In the case 
of the gravitational waves from inspiraling binary black hole systems the effect of the spacetime birefringence is an 
apparent time-dependent change in the inclination angle between the binary system's orbital angular momentum and 
the line-of-sight to the detector. (This change is apparent in the same sense that light is "bent" upon passage nearby a 
strongly gravitating object.) Sufficiently long observations of a binary system will enable this apparent rotation to be 
distinguished from the real rotation caused by spin-orbit and spin-spin angular momentum interactions in the binary 
system. Observations of just this kind will be possible using the LISA gravitational wave detector [3, E3, EBl , which 
will be able to observe the inspiral of massive black hole binaries at redshifts approaching 30 for periods of a year or 
more. Gravitational wave observations of these systems with LISA may thus provide the first test of string theory or 
other quantum theories of gravity: yet another way in which gravitational wave observations can act a unique tool 
for probing the fundamental nature of the universe. 
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